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Abstract 

A measurement is presented of the cross section for D*^ meson production in diffractive 
deep-inelastic scattering for the first time at HERA. The cross section is given for the 
process ep eXY, where the system X contains at least one D*^ meson and is separated 
by a large rapidity gap from a low mass proton remnant system Y. The cross section 
is presented in the diffractive deep-inelastic region defined by 2 < < 100 GeV^, 
0.05 < y < 0.7, xjp < 0.04, My < 1.6 GeV and |t| < 1 GeV^. The D*± mesons are 
restricted to the range pt,d* > 2 GeV and \r]£)* \ < 1.5. The cross section is found to be 
246 lb 54 lb 56 pb and forms about 6% of the con^esponding inclusive D*^ cross section. 
The cross section is presented as a function of various kinematic variables, including 
which is an estimate of the fraction of the momentum of the diffractive exchange carried 
by the parton entering the hard-subprocess. The data show a large component of the cross 
section at low z^'^ where the contribution of the Boson-Gluon-Fusion process is expected 
to dominate. The data are compared with several QCD-based calculations. 
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1 Introduction 



The observation of events with a large rapidity gap in the distribution of the final state hadrons 
at HERA ^ allows the nature of colour singlet exchange in strong interactions to be investi- 
gated. Colour singlet exchange interactions have been successfully modelled ^ in terms of 
phenomenological Regge theory [Q] and, at high energy, are attributed to diffractive or pomeron 
exchange. HERA allows the partonic nature of diffraction to be investigated in deep-inelastic 
scattering (DIS) using the virtual photon as a probe. 

The inclusive diffractive DIS structure function is directly sensitive to the quark content 
of the diffractive exchange [0,^. Information about the gluon content can be inferred indirectly 
from scaling violations. However, the measurement of the hadronic final state in diffraction 
gives further, more direct, information about the gluon content [§|-£p . The production of open 
charm is expected to be particularly sensitive to the gluon content because studies in inclusive 
DIS reveal that the dominant contribution comes from the boson-gluon-fusion (BGF) mecha- 
nism [1T0|]. The presence of the hard scale, provided by the charm quark mass, allows a variety 
of perturbative QCD-based models of diffraction to be tested. 

This article describes the measurement of diffractive open charm production in DIS at 
HERA, which was performed using the HI detector. Measurements of the total D*"^ cross 
section and of differential distributions which explore the dynamics of diffractive charm pro- 
duction are presented. The ratio of the diffractive D*"^ cross section to the inclusive D*"^ cross 
section is also measured. 

The paper is organised as follows. The kinematics of diffractive DIS are introduced in 
section 0. The different theoretical approaches to diffractive charm production are summarised 
in section 0. In section ^ the HI detector, the data selection, the cross section measurement 
procedure and the evaluation of the systematic uncertainties are explained. The results, in the 
form of the total and differential cross sections are presented and discussed in section |[ 



2 Kinematics 

The process studied in this paper is ep — > eXY e(D*'^X')Y and is shown in figure [l]. The 
electron produces a virtual photon 7* (with four-momentum q) which interacts with the proton 
(with four-momentum P). If the interaction takes place via colour singlet exchange, the photon 
and proton dissociate into two distinct hadronic systems X and Y, with invariant masses 
and My, respectively. The system Y is that which is closest to the outgoing proton direction. In 
the case where and My are small compared with the photon-proton centre of mass energy 
W, the two systems are separated by a large rapidity gap. In addition to the standard DIS 
kinematic variables Q"^, y and Bjorken x the following variables are defined 

q.P ' t=^P-Vy)-^ ^= 2q.{P-py) =-^^ 

where py is the four-momentum of Y. The quantity xp may be interpreted as the longitudinal 
momentum fraction, with respect to the incoming proton, of the colourless exchange and t is 
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Figure 1: The process under study in this article is ep eXY e{D*^X')Y. The electron 
(e) couples to a photon ( •y*) which interacts with the proton (p) via colour singlet exchange, pro- 
ducing two distinct final state hadronic systems X and Y. The systems X and Y are separated 
by the largest gap in rapidity in the final state hadrons. 

the squared four-momentum transferred at the proton vertex. In the analysis presented in this 
paper, t and are constrained to be small by the experimental selection and are integrated 
over implicitly. 

In a QCD interpretation in which a partonic structure is ascribed to the colourless exchange 
the lowest order (i.e. 0{a^) ) contribution to the diffractive cross section in the proton infinite 
momentum frame is quark scattering {'j*q ^ q). In this case [3 can be interpreted as the 
fractional longitudinal momentum of the exchange carried by the struck quark. The 0{as) 
contributions are the BGF {'-)*g qq) and QCD-Compton (7*g qg) processes. In the 0{as) 
case the invariant mass squared s of the partons emerging from the hard subprocess is non-zero. 
Therefore, the quantity zjp is introduced 

.r = ^-(n-^) (2) 

which corresponds to the longitudinal momentum fraction of the colourless exchange carried 
by the parton (quark or gluon) which enters the hard interaction. 



3 Models of Diffractive D*^ Production 

A detailed description of the models used in this paper is given in [^. A brief summary focusing 
on the production of open charm in each of the models is given here. For ease of comparison 
with the data in this paper the models are divided into three groups: the 'resolved pomeron' 
model, '2-gluon exchange' models and 'soft colour neutralisation' models. 

In the 'resolved pomeron' model [ ]TT1 ] the diffractive cross section factorises into deep- 
inelastic scattering from the pomeron and a pomeron flux factor, motivated by Regge theory. 
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which describes the probability for finding a pomeron in the proton. Within this picture, the 
partonic content of the pomeron has been determined by QCD analyses of HERA diffractive 
data [@,[1^. The parton distributions, obtained from fits to the data, contain a dominant gluon 



distribution. Open charm is produced in the resolved pomeron model by the BGF process, 
where the photon interacts with a gluon of the pomeron carrying a fraction zjp of the pomeron 
longitudinal momentum. 

In '2-gluon exchange' models diffractive DIS is studied in the proton rest frame by consid- 
ering qq and qqg photon fluctuations as colour dipoles scattering off the proton. Open charm 
can be produced when the photon fluctuates into cc or ccg states. The simplest realisation of 
net colour singlet exchange between these partonic fluctuations and the proton at the parton 



level is a pair of gluons with opposite colour [13]. In perturbative QCD, the cross section for 



2-gluon exchange is related to the square of the /c^-unintegrated gluon density of the proton 
/cy) [[T^,|l5]l, where kx is the parton transverse momentum relative to the proton direction. 



In the 'saturation' model [ [T^ the calculation of the qqg cross section is made under the as- 
sumption of strong kr ordering of the final state partons, where ki^^ ^ ki^'^'' . In an alternative 



approach [[T^,[T8p (hereafter referred to as 'BJLW') the calculation of the qqg final state includes 
configurations without strong kr ordering. In this model all outgoing partons are required to 
have high fcr and the minimum value for the final state gluon transverse momentum k^g is a 
free parameter which can be tuned to describe the data. The 2-gluon exchange calculations are 
performed under the assumption of low xip (xp < 0.01) to avoid contributions from secondary 
reggeon exchanges which correspond to quark exchange in the models. 

In 'soft colour neutralisation' models an alternative approach to diffractive DIS is given 
which leads to very similar properties of inclusive and diffractive DIS final states. In the soft 



colour interaction model (SCI) [[191], open charm is produced via BGF from the gluon distri- 
bution of the proton. It is then assumed that the partons produced in the hard interactions 
can exchange soft gluons with the background colour field of the incoming proton leaving all 
momenta unchanged. Large rapidity gap events may be produced this way when the soft in- 
teractions lead to a net colour singlet exchange. The probability for soft colour exchange is 
assumed to be independent of the kinematics of the hard scattering process. 



The generalised area law (GAL) [ pOQ approach is a modification of the Lund String model 
[PH]. The production mechanism for open charm is similar to that in the SCI model except 
that it is formulated in terms of interactions between the colour strings connecting the partons 
in an event. In this model the probability for a soft colour interaction is not constant but is 
exponentially suppressed by the difference between the areas in momentum space spanned by 
the strings before and after the colour rearrangement. 



The 'semi-classical' model [ ^zj ] is a non-perturbative model based on the dipole approach. 
In the proton rest frame the photon fluctuations scatter off a superposition of soft colour fields 
representing the proton. In this approach, the qqg fluctuation is expected to be dominant for 
open charm production [O]. If the gluon is the lowest k^ parton, then the contribution can be 



related to BGF in the proton infinite momentum frame. 

Comparison of the data with the 'resolved pomeron', the '2-gluon exchange' and the 'semi- 
classical' models is facilitated by their implementation within the RAPGAP Monte Carlo gen- 



erator [p4[]. The predictions of the SCI and GAL models are calculated using the AROMA 
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Monte Carlo generator [p5]]. The cross section predictions in this article are all calculated as- 
suming a charm quark mass nic = 1.5 GeV. For the hadronisation fraction /(c D*^) the 
value 0.233 ± 0.010 ± 0.011 is used. The momentum fraction of the charm quark carried 
by the D*"^ is calculated using the Peterson et al. model with the fragmentation parameter 
e = 0.078. 



4 Experimental Procedure 

The data presented in this analysis were collected over the years 1996 and 1997, when HERA 
collided positrons with energy = 27.5 GeV with protons of energy Ep = 820 GeV. Re- 
quiring all essential detector components to be operational the available integrated luminosity 
is IQ.lpb"^. Further details of this analysis beyond those given here can be found in [^. 



4.1 The HI Detector 

A short overview of the detector components most relevant for the present analysis is given 
here. A detailed description of the HI detector can be found in [p9|]. The z-axis of the HI 
detector is taken along the beam direction such that positive z values refer to the direction of 
the outgoing proton beam, referred to as the 'forward' direction. 

Charged particles emerging from the interaction region are measured by the central track- 
ing device (CTD) in the range —1.5 < r] < 1.50. The CTD comprises two large cylindri- 
cal central jet drift chambers (CJC) and two z chambers situated concentrically around the 
beam-line within a solenoidal magnetic field of 1.15 T. The resolution achieved by the CTD 
is o-(pj')/pt — O.Olpr/GeV. The CTD also provides triggering information based on track 
segments in the r — (p plane from the CJC and the position of the vertex using a double layer of 
multi-wire proportional chambers (MWPC). The energies of final state particles are measured 
in the Liquid Argon (LAr) calorimeter which surrounds the tracking chambers and covers the 
range —1.5 < r] < 3.4. The backward region (—4.0 < rj < —1.4) is covered by a lead- 
scintillating fibre calorimeter (SPACAL [P^]) with electromagnetic and hadronic sections. In 
front of the SPACAL, the Backward Drift Chamber (BDC) provides track segments of 
charged particles. 

Detectors close to the beam pipe in the direction of the outgoing proton are used in the se- 
lection of large rapidity gap events. These are the Forward Muon Detector (FMD), the Proton 
Remnant Tagger (PRT) and the Plug calorimeter (PLUG). The FMD is located dX z = 6.5 m 
and covers the pseudorapidity range 1.9 < rj < 3.7 directly. The PLUG allows energy mea- 
surements to be made over the range 3.5 < rj < 5.5. Particles produced at larger rj can also 
be detected because of secondary scattering with the beam-pipe. The PRT, a set of scintillators 
surrounding the beam pipe at 2 = 26 m, can tag hadrons in the region 6.0 "^rj < 7.5. 

'The pseudorapidity 7/ of an object detected with polar angle is defined as 77 = — In tan(6'/2). 
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4.2 Event Selection 



The events were triggered by an electromagnetic energy cluster in the SPACAL, in coincidence 
with a charged track signal from both the MWPC and the CJC. The positrons are identified 
in the SPACAL as clusters with energy E'^ > 9 GeV which have properties consistent with 
electromagnetic deposition, and for which the centre of gravity of the cluster matches a charged 
track segment in the BDC to within 2.5 cm. The selected events are also required to have a 
reconstructed vertex from the CTD within ±35 cm of the nominal vertex. In order to suppress 
events with initial state photon radiation the summed E — of the event calculated using all 
reconstructed final state particles, including the positron, is required to be greater than 35 GeV. 
The kinematic region covered by the measurement is 2 < < 100 GeV^ and 0.05 < y < 0.7. 
To minimise the correction due to QED radiative effects, Q^, y and x are reconstructed from the 
energy and angle 9' of the scattered positron and the hadronic final state using the 'S method' 



Diffractive events are selected experimentally by the absence of activity in the outgoing 
proton region. No signal above noise thresholds is allowed in the FMD, the PRT, the PLUG and 
the most forward part (r] > 3.3) of the LAr calorimeter. This ensures that there is a large rapidity 
gap covering at least 3.3 < r] < 7.5 between the photon dissociation system X and the proton 
remnant system Y. Monte Carlo studies show that the absence of particles in the detectors close 
to the beam pipe restricts the mass of the proton remnant system to My < 1.6 GeV and the 
momentum transfer to the proton to |t| < 1 GeV^. 

The four-momentum of the system X, which is well contained in the central detector, is 
reconstructed using information from the LAr and SPACAL calorimeters together with the CJC 



[33]. The variable xjp is calculated from 



X. = + (3) 

2Ep ^ 

where E and are the energy and longitudinal momentum of each final state particle in the 
laboratory frame, and the sum runs over the scattered positron e' and all detected particles in the 
photon dissociation system X. The quantity /3 is calculated from /3 = x/xjp. The cross section 
is restricted to the range xjp < 0.04 to suppress contributions from non-diffractive scattering 
and secondary reggeon exchanges. 

In this paper an hadronic observable z^^ is constructed which is analogous to x°^'^ for in- 
clusive D*"^ production which was measured in [Q]. In the resolved pomeron picture is an 



approximation to the momentum fraction 2;jp of the pomeron carried by the interacting gluon 
(see equation H). The observable is defined as 



(4) 



where M^g is a hadron level estimate of s which is constructed from the scattered positron and 
the D*^ meson in an identical manner to that used for the gluon momentum fraction x°g^ in Q] . 
Monte Carlo simulations show that the resolution in the hadronic variable Zp'^ is approximately 
30%, and that there is a good correlation between Zp^ and zp as calculated from the kinematics 
of the outgoing partons. The variable z^^ can be interpreted as the fraction of the energy of the 
system X which is carried by the cc pair emerging from the hard scattering. 
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4.3 Reconstruction of D*^ Mesons 



The D*^ mesons are reconstructed using the D*^ 
decay channel 



D mass difference method [35] in the 



D 



*± 



7r 



slow 



7r 



slow 



+C.C. 



(5) 



which has a branching fraction of 2.59% The reconstruction method is detailed in [[37|]. 
The decay products are detected in the CTD and are required to have a transverse momentum 
Pt of at least 140 MeV for the tTsIow and 250 MeV for both the K and vr. 




Figure 2: Distribution of the mass difference AM = M(K^n'^nfi^^) — M^K^n"^), with a 
curve fitted to the form a{AM — M^)''+ Gaussian. 

The invariant mass of the Ktt combination has to be consistent with the -D° mass within 
±80 MeV. After cuts on the direction (|?7(_ft'7r7r) I < 1.5) and transverse momentum (^^(KTrTr) > 
2 GeV), the mass difference distribution AM = M{KtxtXsIow) — M^Ktt) is plotted in figure |^. 
The number of D*"^ candidates is determined by fitting the histogram in figure ^ with a Gaussian 
distribution for the signal plus a background function a{AM — M^)^, where M^ denotes the 
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mass of the pion. The position and width of the Gaussian are fixed to values taken from a higher 
statistics sample of events where no diffractive cuts were applied The normalisation of the 
Gaussian and the background parameters a and b are allowed to vary. The resulting number of 
detected D*^ mesons is 46 ± 10. 



4.4 Cross Section Measurement 

Monte Carlo simulations are used to correct the data for the effects of losses and migrations 
due to the finite resolution of the HI detector. The efficiency is calculated by running the 
HI detector simulation program on a sample of D*^ events from the diffractive Monte Carlo 
generator RAPGAP in the resolved pomeron mode with the t dependence of the cross 
section parameterised as e^^'*'. The RAPGAP program is used to model events which contain 
an elastic proton (My = rup) in the kinematic range xp < 0.1. Migrations from xp > 0.1 
or from large values of My (My > 5 GeV) are modelled by using a simulation of the heavy 



quark generator AROMA [ PSQ in the inclusive mode. The contribution is of the order 5% to 
the selected sample of events. An additional correction of —8% ± 6% is applied to account for 
the net smearing across the = 1.6 GeV boundary. Since only elastically scattered protons 
have been simulated in RAPGAP, this correction is evaluated using the proton dissociation 
simulation in the DIFFVM [ p8| ] generatoiQ. A further correction of +4%±1% takes into account 
diffractive events rejected due to fluctuations in the noise level in the FMD. This correction is 
estimated directly from the data, using a sample of randomly selected events, not correlated 
with a physics trigger. An additional source of background is the contribution of reflections in 
the mass window, coming from channels other than that defined in equation ^. They are 



estimated, from simulations using the AROMA Monte Carlo, to be 3.5% Pq]. The contribution 



from photoproduction background is found to be negligible. QED radiative corrections were 



calculated to be approximately 2% using the RAPGAP program interfaced to HERACLES \39]. 



4.5 Systematic Uncertainties 

The following sources of systematic error are taken into account 

• The uncertainty in the physics model for D*^ production used to compute the efficiency 
corrections is estimated by varying the shapes of the kinematic distributions in the sim- 
ulations beyond the limits imposed by previous measurements or the present data. This 
is done by reweighting the xp distribution to that observed in data; the P distribution by 
(1 ± 0.3/3) and the t distribution by e^"^*. The resulting systematic uncertainties on the 
cross section measurements range between 10% and 20% with the largest contribution 
originating from the variation of the xp distribution. The uncertainties are verified using 
simulations of models with different underlying kinematic distributions. 

• The total uncertainty due to the reconstruction efficiency, mass and momentum resolution 
of the central tracker for the three tracks was estimated in the analysis of the inclusive DIS 



D*"^ cross section to be '^l^ [BTI]. 



^For the correction, it is assumed that the ratio of diffractive proton elastic to diffractive proton dissociative 
interactions is 1 : 1. 
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An error of 8% is found by varying the details of the fitting procedure used to obtain the 
number of D*^ mesons. 

The uncertainty in the correction due to the smearing of events across the boundary My = 
1.6 GeV is estimated by varying in the DIFFVM simulation: the efficiency of the forward 
detectors, the assumed My distribution, the ratio of double to single dissociation between 
0.5 and 2 and the assumed t dependence for double dissociation. This contributes 6% to 
the systematic error. 

The uncertainty in the trigger efficiency gives a contribution of 5% to the systematic error. 

The uncertainty due to the assumed charm fragmentation scheme is estimated by using 
parameterisations of the Peterson model and the standard Lund string model in JETSET 



[plQ. This leads to an average uncertainty of 5% in the cross sections. 



The uncertainty in the correction due to QED radiative effects is estimated as 3%. 

The number of events migrating into the sample from xp > 0.1 or My > 5 GeV is varied 
by ±50%, leading to an average systematic error of 3%. 

Other sources of systematic error are the uncertainty in the measured energy and angle of 
the scattered positron, uncertainties in the hadronic energy scale of the liquid argon and 
SPACAL calorimeters, the uncertainty in the luminosity measurement, the uncertainty on 
the fraction of events lost due to noise in the FMD and the uncertainty in the branching 
ratio for the measured decay channel. Each of them is responsible for an error of no more 
than 2.5%. 



The total systematic error for each point has been obtained by adding all individual contri- 
butions in quadrature. It ranges between 20% and 30% and for most data points is similar in 
magnitude to the statistical error. 



5 Results 

The total D*"^ production cross section for the kinematic region 2 < < 100 GeV^, 0.05 < 
y < 0.7, X]p < 0.04, My < 1.6 GeV, |t| < 1 GeV^, pt,d* > 2 GeV and |r/B. | < 1.5 is 

a{ep e{D*^X')Y) = 246 ± 54 (stat.) ± 56 (syst.) pb . (6) 

The ratio of the diffractive D*^ cross section to the inclusive D*"^ cross section measured 
in the same kinematic range defined in terms of Q"^, y, Pt,d* and r]D* is found to be 

5.9 ± 1.1 (stat.) ± 1.1 (syst.)%, (7) 

where the inclusive D*^ cross section has been determined as in [P7|]. The error in the ratio is 
dominated by the uncertainties pertaining to the measurement of the diffractive cross section. 
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Model 


cross section (pb) 


resolved 
pomeron 


HI fit 2 
HI fits 
ACTW fit D 


368 
433 
481 


soft colour 
neutralisation 


SCI 

GAL 
semi-classical 


203 
328 
196 


HI ] 


Data 


246 ± 54 (stat.) ± 56 (syst.) pb 



Table 1: The predictions for the total dijfractive D*^ cross section for two groups of models: 
the resolved pomeron and soft colour neutralisation approaches. The bottom row shows the 
cross section measured in the data. 

In table [I] the total cross section is compared with some of the phenomenological models 
discussed in section ^. The first three rows of table [1] show the predictions for the cross section 
for three different sets of parton parameterisations within the resolved pomeron model. The first 
two predictions are based on the parton distributions of the pomeron and sub-leading exchange 
from the leading order DGLAP analysis of F2 from HI The 'ACTW fit D' parameterisation 
is the best combined fit in [[T|] to HI and ZEUS data and ZEUS diffractive dijet data. All 
three sets of parton parameterisations give acceptable descriptions of F2 . AH of the predictions 
using the three parton parameterisations exceed the data although the parameterisation with the 
flat gluon distribution ('HI fit 2') is closest. The predictions shown are calculated with the 
factorisation and renormalisation scales set to /i^ = /Xj = /x^ = + + 4m^. Changing 
this scale to p\ + 4m^ produces an increase of around 20% in the predicted cross sections. 
Similarly, the variation of the charm quark mass by ±0.1 GeV leads to an uncertainty of =f10% 
in the cross sections. Changing e in the Peterson model from 0.078 to 0.035 and to 0.1 produces 
an uncertainty in the cross section predictions of ^5^%. The values shown in the table are 
calculated with Aqcd = 0.20 GeV and the number of active quark flavours in the first order 
expression for is Nj = 4. Selecting Aqcd = 0.25 GeV and Nf = 5 leads to an increase of 
about 10% in the cross sections. The contribution of D*^ production from meson exchange in 
the predictions is less than 7%. 

The cross section predictions from the semi-classical, SCI and GAL models are also shown 
in table [I] and are in agreement with the data. However, none of these models can simultaneously 
reproduce the shapes and normalisations of the differential dijet cross sections [|^]. The semi- 
classical model prediction was calculated using the same factorisation scale as for the resolved 
pomeron model. The SCI and GAL model predictions use /x^ = + 2p|, + 2m^. For each 
of the three models, the uncertainty in the predictions due to the variation of the factorisation 
scale, nic and e are similar to those for the resolved pomeron model. 

Differential cross sections are shown in figures ^ and ^. They represent average values over 
the intervals shown in the figures. 

In figure |3| the cross section is shown as a function of xjp, log^g P ^^'^ ^1p^ compared to 
the QCD-based models described in section]^. The cross sections differential in xp and logio P 
are flat within experimental errors. Figure |3| shows that about 60% of charm production is in 
the region Zp'^ < 0.2 where the contribution of the BGF process is expected to dominate. 
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Figure 3: Cross sections a{ep — > e{D*'^X')Y) as a function of (a) xp; (b) logiQ /? and(c) Zp^. 
The data are points with error bars ( inner: statistical, outer: total). Each distribution is plotted 
three times to allow comparison with the three groups of models described in the text. 
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Figure 4: Cross sections a{ep — > e{D*^X')Y) as a function of (a) logioQ^; (b) p^^o* (^) 
r)D*- The data are points with error bars (inner: statistical, outer: total). Each distribution is 
plotted twice to allow comparison with two of the groups of models described in the text. 
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The discrepancy between the resolved pomeron model predictions and the data is pro- 
nounced in the low xp, high P and high z^^ regions which are all correlated with low values of 
M^. Of the three parton parameterisations the calculations which use 'HI Fit 2' are shown to 
come closest to the data in this region. All three parameterisations are consistent with the data 
in the high region. 

The two gluon exchange models, directly applicable in the low xp region (xjp < 0.01) 
are also compared to the data in figure |3|. The sensitivity of this measurement to the value of 
the transverse momentum cut-off kj^g of the gluon in the qqg state within the BJLW model is 
also studied. Calculations which use cut-off values of 1.0 GeV and 1.5 GeV both give a fair 
description of the data in the low region (i.e. the low xjp, high P, high 2;^* domain). These 
data also offer sensitivity to the relative contribution of the scattering of the qq fluctuation, which 
is shown as a shaded zone, and forms a sizeable component of the total two gluon exchange 
cross section. The saturation model reproduces reasonably well the normalisation of the data in 
the low xjp range, in which it is expected to be applicable, but also provides a good description 
of the data in the remaining region of phase space. 

The semi-classical model gives a good description of the distributions shown in figure |^. 
Both the SCI and GAL models provide a satisfactory description of the spectra although the 
GAL model tends to overestimate the data in the low domain. 

In figure ^ the D*^ cross section is plotted differentially at 2 values of logio Q^^ Pt d* 
r]D* , where j^. is the transverse momentum of the D*^ in the 7*p centre of mass system. The 
data tend to fall off with higher values of each of these variables. 

Since these distributions are integrated over the full xjp range of this measurement, a com- 
parison is made only with resolved pomeron calculations and soft colour neutralisation models. 
Both sets of models provide a reasonable description of the data. 

6 Conclusion 

The dynamics of open diffractive charm production in DIS have been studied for the first time 
at HERA. The total D*"^ production cross section in the kinematic range 2 < < 100 GeV^, 
0.05 <y < 0.7, xjp < 0.04, My < 1.6 GeV, \t\ < 1 GeV^, pt,d* > 2 GeV and \r]D'-\ < 1.5 
has been found to be 246 ± 54 (stat.) ± 56 (syst.) pb. In the studied region about 6% of the 
total D*"^ cross section is produced diffractively. 

The cross section has been measured as a function of xp, log^Q (3, Zp^, log^g Q^^ Pt,d* 
r]D*. The data show a sizeable component of charm production in the low Zp^ region which is 
suggestive of the dominance of the contribution from the boson-gluon-fusion process. 

A number of QCD-based models which give a good description of the inclusive diffractive 
cross section were compared with the measurement. A reasonable description of the data is 
provided by a model based on the resolved pomeron picture using various assumptions for the 
partonic composition of the colourless exchange. A parton parameterisation containing a flat 
gluon dependence ('HI Fit 2') comes closest to the data. Predictions of two gluon exchange 
processes were found to match the data in the low xp region. Soft colour neutralisation models 
give a satisfactory description of the data. 
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